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We report about two new rigorous results on the non-analytic properties of thermodynamic po- 
tentials at first order phase transition. The first one is valid for lattice models (d > 2) with arbitrary 
finite state space, and finite-range interactions which have two ground states. Under the only as- 
sumption that the Peierls Condition is satisfied for the ground states and that the temperature 
is sufficiently low, we prove that the pressure has no analytic continuation at the first order phase 
transition point. The second result concerns Ising spins with Kac potentials J 1 {x) = j d ip(-yx), where 
< 7 < f is a small scaling parameter, and ip a fixed finite range potential. In this framework, 
we relate the non-analytic behaviour of the pressure at the transition point to the range of inter- 
action, which equals 7 _1 . Our analysis exhibits a crossover between the non-analytic behaviour of 
finite range models (7 > 0) and analyticity in the mean field limit (7 \ 0). In general, the basic 
mechanism responsible for the appearance of a singularity blocking the analytic continuation is that 
arbitrarily large droplets of the other phase become stable at the transition point. 
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INTRODUCTION 
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The first theory of condensation originated with the 
celebrated equation of state of van der Waals 0] : 



P 



RT. 



(1) 



When complemented with the Maxwell Construction 
(or "equal area rule"), Q leads to isotherms describing 
general characteristics of the liquid-vapor equilibrium, 
including the existence of a critical temperature. The 
isotherms obtained with the van der Waals-Maxwell 
Theory have a very simple analytic structure: they are 
analytic in a pure phase and have analytic continua- 
tions along the liquid and gas branches, through the 
transition points. These analytic continuations, which 
were originally interpreted as describing the pressure of 
metastable states, are provided by the original isotherm 
given in 

The theoretical question of knowing whether the 
results predicted by the van der Waals Theory can be 
derived from first principles of Statistical Mechanics 
remained a longstanding problem during a large part 
of the twentieth century. The theories of Mayer [2j and 
Yang-Lee were decisive contributions to the theory of 
phase transitions, but didn't give an answer concerning 
the delicate question of the analytic continuation at 
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transition points. With regard to this latter property, 
two scenarios were discussed in the fifties and sixties. 

The first one was essentially based on the mean 
field (or Bragg- Williams) approximation £|. In this 
approach, the interaction is replaced by an infinite range 
and infinitely weak potential. The central characteristic 
of the effective model obtained after this approximation 
is that the spatial positions of the particles don't play 
any role. As a consequence, an exact computation of 
the partition function leads to the same behaviour as in 
the van der Waals-Maxwell Theory: at low temperature, 
thermodynamic potentials are analytic in a pure phase, 
and have analytic continuation at transition points. 
Katsura [f| conjectured that this scenario holds also for 
short range models, like the Ising model (see also the 
discussion below). 

The second argument, totally different in spirit, 
originated with the so called "droplet mechanism" of 
the condensation phenomenon, proposed by Andreev 
@, Fisher Q and Langer ||. This mechanism, as 
opposed to the mean field approximation, predicts that 
the finitcness of the range of interaction plays a crucial 
role in the analytic properties of the thermodynamic 
potentials. Namely, when the range of interaction is 
finite, droplets of any size are stable at the condensation 
point, and although the probability of occurence of large 
droplets is very small, they yield a contribution of the 
order k\~z= T to the fc-th derivative of the pressure. Kunz 
and Souillard were led to the same conclusions after 
having studied a similar model, related to percolation ||. 
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Subsequent papers on the subject, in which no definite 
answer was given, include luJ, More recent 

studies can also be found in 



RIGOROUS RESULTS 

The first rigorous result was the study of Isakov |l6j on 
the Ising model, which confirmed the droplet predictions: 

Isakov (1984): In dimension d > 2, at low enough 
temperature, the pressure of the Ising model in a mag- 
netic field A, p = p{\), is infinitely differ entiable at 
X = 0^, but has no analytic continuation from {A < 0} 
to {A > 0} accross A = 0, or vice versa. 

This result is obtained by showing that the derivatives 
of the pressure at A = behave like 



P 



(0 ± ) ~ C k k\^ . 



(2) 



In a second paper , Isakov tried to extend this result 
to general two phase lattice models. He had, however, 
to introduce hypothesis that are not easy to verify in 
concrete models. We now present our results. 



Let < 7 < 1 be a small scaling parameter, and consider 
the Kac potential J 7 (x) = 7^(72:), together with the 
hamiltonian 



H 



(5) 



Let / 7 = / 7 (m) denote the free energy of this model, 
with fixed magnetization m G [—1, +1]. The Theorem of 
Lebowitz-Penrose [2lJ gives a closed form to the free en- 
ergy in the van der Waals limit 7 \ (called sometimes 
the Kac or mean field limit), and justifies the Maxwell 



construction. Let fo(m) 
ureQJ 



lim 7 \ N o f~f(m). Then (see Fig- 



fo(m) 
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nvclope of j — —m 2 — — J(m)| , (6) 



where I(m) equals 



r , , 1 — m, 1 — m Urn, 1 + m ._. 
Hm) = ^— log — = 7T~ l °S — s~ • ( 7 ) 



Two Phase Models. Consider a lattice model with fi- 
nite state space at each site of Z d , d > 2. Let H be a 
hamiltonian with finite range periodic interaction, hav- 
ing two periodic ground states ipi , ip% , so that the Peierls 
Condition is satisfied ^^|. Let V be a periodic poten- 
tial with finite range interaction, so that the perturbed 
hamiltonian 



H\ — Hn 



XV 



(3) 



splits the degeneracy of Hq. That is, H\ has a single 
ground state ip2 when A < and a single ground state 
-01 when A > 0. Denote by p = p(X) the pressure of 
the model (at inverse temperature (3). Let 8 > 0. The 
general theory of Pirogov-Sinai guarantees that if j3 
is large enough, then there exists A*(/3) £ (—6, +6) such 
that the pressure has a first order phase transition at 
A*(/3). Our first result [2(| is the following: 

Theorem 1 There exists /3q > such that for all 
P > Po> the pressure is analytic in A on (—5, X*(f3)) 
and (X*(/3),+8), but has no analytic continuation from 
(—5, X*(f3)) to (X*(f3),+5) across X*(f3) or vice-versa. 

Kac Potentials and the van der Waals Limit. Con- 
sider an Ising fcrromagnet, with a spin <Tj G {+1,-1} at 
each site of Z d , d > 2. Let ip : R d -> M+, supported by 
the cube [— such that 



ip(x)dx = 1 . 



(4) 



FIG. 1: The free energy in the van der Waals Limit. 

When f3 > 1, fo(m) has a plateau [— m*{[3), +m*(/3)], 
where m*{[3) is the positive solution of the mean field 
equation m = tanh(/3m). As a consequence of the 
Lebowitz-Penrose Theorem, all the analytic properties 
of the free energy are known explicitly after the van der 
Waals limit: fo is analytic on the branches (—1, — m*(/3)) 
and (+m*(/3), +1), and has analytic continuation along 
the paths m / —m*([3), m \ +m*(/3). The analytic 
continuation, which is unique, is given by the mean field 
free energy —\rn 2 — ^7(m). After the van der Waals 
limit, the scenario is thus the same as in the van der 
Waals-Maxwell Theory. 

Consider the specific choice f{x) = 2~ d l(x), where l(-) 
is the indicator of the cube: l(x) = 1 if x G [— 
otherwise. For a fixed < 7 < 1, J 7 is finite range 
and Theorem ^ can be used, but only for temperatures 
P > ^0(7), with lhn 7 \o/3o(7) = +00. Our result j2^| is 
given hereafter. It holds at low temperature, uniformly 
in the range of interaction. 

Theorem 2 There exists (3q, independent 0/7, such that 
for all (3 > Pq and all < 7 < 1, the free energy / 7 
is analytic on (— 1, — m*{f3, 7)) and (+m*(f3, 7), +1), but 
has no analytic continuation along the real paths m /* 
-m*(/3,7), m \ +m*(f3,j). 

This result shows that, as opposed to the mean field 
behaviour, finite range interactions, even of very long 
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range, imply absence of analytic continuation at tran- 
sition points. A crucial ingredient for the proof of 
Theorem [21 is the use of a coarse-graining technique 
proposed by Bovier and Zahradmk it allows to 

obtain a range of temperature that is uniform in 7. 

We study the pressure p 7 = p 7 (A), in which the con- 
straint on the magnetization is replaced by a magnetic 
field A. The pressure and free energy are related by a 
Legendre transform: 

/ 7 (m) = sup (hm — p 1 (A)) . (8) 

A 

By the Theorem of Yang-Lee, p 7 is analytic in A on 
{A < 0} and {A > 0}. Our main result is a precise char- 
acterization of the properties of p 7 along the path A \ 
(using symmetry, we need only consider fields A > 0). 

Theorem 3 There exists [3q, independent ofy, such that 
for all (3 > Po an d a ^ 1 > 0, all the limits p^ (0 + ) = 
limx\oP7 (A) exist, but pressure has no analytic contin- 
uation from {A > 0} to {A < 0} accross A = 0. More pre- 
cisely, there exists integers ki(j), ^2(7), ^1(7) < ^2(7); 
with lim 7 \^o ^1(7) = +00, such that 

|p 7 fe) (0 + )l < Cffc! when k < , (9) 

|pW( +)| > C|jfc!^r w hen k > fc 2 ( 7 ) . (10) 

The constant C\ is independent of 7 and k, C2 = 
C 2 ( 7 ,/3)>0, and fc x ( 7 ) =7^. 

That is, the large order derivatives reveal the non- 
analytic feature of the singularity, although a signature 
of the mean field (analytic) behaviour can be detected 
in the low order derivatives. We have illustrated this 
crossover on Figure [21 

p 7 fc) (0±)~fe! ^'(OVfclA 



fci(7) k 2 (j) 
FIG. 2: The crossover in the derivatives of the pressure. 



METHOD 

The pressure has a singularity only in the thermo- 
dynamic limit. However, we study the system in large 
finite volumes, and obtain bounds on the derivatives 
that are uniform in the volume. At the end we prove 
that it is possible to interchange the operations of taking 
the derivative and the thermodynamic limit. 

The method used to obtain lower bounds on the deriva- 
tives of the pressure at finite volume is inspired by the 



technique of Isakov. Let A be a finite cube in Z d with a 
fixed boundary condition, and Z(A) be the correspond- 
ing partition function. One enumerates all possible con- 
tours inside A: ri, T2, ■ ■ ■ , T n , in such a way that 
V( r ») ^ v ( r j) wnen i < j (V(Ti) denotes the volume 
of the interior of the contour Ti). One then defines the 
restricted partition functions Zi(A), i = 0,...,n. By 
definition, Zq (A) is the partition function computed for a 
system containing no contours, and Zi(A) is the partition 
function computed for a system containing no contour Tj 
with j > i. Obviously, 

Z(A) = Z (A)f[^j- (11) 

For the proof of Theorem^ there is only the ground state 
configuration contributing to Zq(A). For the proof of 
Theorem[21 Zq (A) is the partition function of a restricted 
phase, describing small local fluctuations of the ground 
state. Let 

Notice that we have the fundamental relation 

Z i (A) = 2i_i(A) + 2?_ 1 (A), (13) 

where the contour Tj appears in each configuration con- 
tributing to Z*_ 1 (A). A precise analysis of the phase 
diagram shows that A 1— ► ma(Fj)(A) is analytic in a disc 
Ui centered at A = A*(/3) (resp. A = for the Kac ferro- 
magnet), with a radius of order V(Ti)^' 3 . In the domain 
Ui, u^iTi) can be represented as follows: 

u A (ri) = log (l + |M^) = Iog(l + e 9A(r * } ) ■ (14) 

The dependence of g\(Ti) on the volume A is weak. 
Moreover, g\(Ti) can be decomposed into a surface term 
and a volume term, like in the droplet model. Then, by 
choosing a path of integration C C Ui, 

^(A) U (A) = ^-y c(A „ A , )fc+1 dA. (15) 

The observation of Isakov is that u\(Ti) ~ e 9A(r\) Qn 
Ui, and that for a given large enough k, one can es- 
timate precisely the Cauchy integral (|15|) . for all large 
enough contours, by a stationary phase method, choos- 
ing suitably the path of integration C. In this way one 
gets a contribution to the k-th derivative of the pressure 
of order A k kl~ s ^ rj . For the other contours, only an upper 
bound can be obtained on the integral, of the same or- 
der B k k\' 3 ^ T . The crucial point is therefore to have large 
enough neighbourhoods Ui, in order to show that A > B. 
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DISCUSSION 

In the framework of Kac potentials, the role played 
by the range of interaction in the analyticity properties 
of the pressure can be clarified by the following discus- 
sion. When A > 0, our analysis allows to decompose the 
pressure in two distinct parts: 



P^f — ^* 7 ~\~ Q-y 



(16) 



On one hand, r 7 is constructed with the partition func- 
tion Zq(A) of (11) , and describes a homogeneous phase 
with positive magnetization, containing no droplets of 
the — phase. When 7 \ 0, r 7 converges to the pressure 
of the mean field model. On the other hand, g 7 contains 
the contributions from the droplets of the — phase, which 



are all stable at A = 0, and g 7 = 0( 



-P-y- 



Namely, 



the main contribution to q 7 comes from the smallest 
droplets, which live on a coarse-grained lattice whose 
cells have side length 7 -1 . Then, the pressure r 7 behaves 
analytically at A = 0, i.e. ^^(O^) ~ fc! for all k, but 
q~f is responsible for the absence of analytic continuation 

at A = 0, since q- y ( - k \0 ± ) ~ fc! 3 ^ 1 for large enough 
k. The combination of these two behaviours leads to 
a crossover in the derivatives, as was shown in Theorem[3] 

Our results also have an important consequence re- 
arding the theory of condensation initiated by Mayer 
In this theory, the pressure of a non-ideal gas is de- 
scribed, near z = (z is the fugacity), by a convergent 
Taylor expansion, given by the Mayer series: 



iz 



(17) 



i>i 



(The bi are the cluster coefficients.) That is, the pressure 
is analytic in the gas phase. Then, the condensation 
point is defined to be the first singularity encountered 
when Ijl7|) is continued analytically along the positive 
real line z > 0. It was suggested by some authors 
that this method could actually lead to a wrong 
determination of the condensation point: the analytic 
continuation of the Mayer series might not "see" the 
real transition point z*, situated somewhere between 
and z* M : < z* < z* M . This is indeed the case in the 
mean field approximation: the system docs not "see" 
the condensation point, since there are no droplets of 
the liquid phase inside the gaseous phase. We saw 
that if one suppresses the condensation mechanism by 
retaining, in - only the term r 7 , then there is an 
analytic continuation of the pressure. Our analysis 
shows that the phenomenon of condensation itself 
is responsible, through the singular term g 7 , for the 
singularity. Therefore, the latter method of determining 
the condensation point gives the correct result, at least 
for a large class of lattice gas models. 



To conclude, we mention that the problem of know- 
ing whether the pressure can be continued analytically 
around the singularity, in the complex plane, remains 
open. It is not clear, in this case, whether the droplet 
models can be used as a guiding mechanism, even to give 
a heuristic description [25| ■ 
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